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Abstract. Recently, in [T| we developed a parametric reconstruction method to 
a homogeneous, isotropic and spatially flat Friedmann-Robertson- Walker (FRW) 
cosmological model filled of a fluid of dark energy (DE) with constant equation of 
state (EOS) parameter interacting with dark matter (DM). The reconstruction method 
is based on expansions of the general interaction term and the relevant cosmological 
variables in terms of Chebyshev polynomials which form a complete set orthonormal 
functions. This interaction term describes an exchange of energy flow between the DE 
and DM within dark sector. In this article, we reconstruct the interaction function 
expanding it in terms of only the first four Chebyshev polynomials and obtain the 
best estimation for the coefficients of the expansion assuming three models: (a) a DE 
equation of the state parameter w = —1 (an interacting cosmological A), (b) a DE 
equation of the state parameter w — constant with a dark matter density parameter 
fixed, (c) a DE equation of the state parameter w = constant with a free constant dark 
matter density parameter to be estimated, and using the recent five test data: Union2 
SNe la data set from "The Supernova Cosmology Project" (SCP) composed by 557 
type la supernovae, the Baryonic Acoustic Oscillations (BAO), the CMB anisotropics 
from 7-year WMAP experiment, the Hubble expansion rate data and the X-ray gas 
mass fraction. In all the cases, the preliminary reconstruction shows that in the best 
scenario there exist the possibility of a crossing of the noninteracting line Q = in 
the recent past within the Icr and 2a errors from positive values at early times to 
negative values at late times. This means that, in this reconstruction, there is an 
energy transfer from DE to DM at early times and an energy transfer from DM to 
DE at late times. We conclude that this fact is an indication of the possible existence 
of a crossing behavior in a general interaction coupling between dark components. 
Finally, we conclude that in this scenario, the observations put strong constraints on 
the strength of the interaction so that its magnitude can not solve the coincidence 
problem or at least alleviate significantly. 
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1. Introduction 

Recent observations of tlie apparent magnitude of Supernova la (SNIa) suggest that 
the universe is in a stage of recent acceleration [2]- [3]. This has now been confirmed 
by another sets of independent observational data as such the measurements of the 
galaxy power spectrum and the Baryon Acoustic Oscillations (BAO) detected in the 
large-scale correlation function of luminous red galaxies in the experiment Sloan Digital 
Sky Survey (SDSS) |1]-[H], the power spectrum of the cosmic microwave background 
(CMB) anisotropies measured in the experiment Wilkinson Microwave Anisotropy 
Probe (WMAP) [9]-[TQ]. The most simplest explanation to the recent acceleration 
is the assumption of the existence of a cosmological constant A assumed to stand 
for the vacuum density energy (pA ^ A/8ttG). However, its inferred value from 
observations p'^^ ~ (10~^^Gev)^ is too tiny when is compared with the theoretical value 
= (lO^^Gev)^ estimated from the application of the quantum field theory to the 
Planck scale. This is known as the cosmological constant problem [II]-[33]. In addition 
to this problem, the cosmological constant model has a new one known in the literature 
as the cosmic coincidence problem [13], [31]-|40] which can be established as: why are 
both dark densities (dark matter and dark energy) of the same order of magnitude at the 
present whereas that they were so different in most of the past evolution of the universe?. 
In order to avoid the problem of fine tuning of the cosmological constant, sometimes 
it is assumed a null contribution of the vacuum energy density to the cosmological 
constant which is fixed to zero and therefore it does not contribute to the gravitational 
sector. In this scenario, we need a different explanation for the present acceleration of 
the universe. In this direction, it has been proposed a new form of matter-energy named 
dark energy (DE) which can be a perfect fluid [3], [8], [10], |ll]-|l2] or a slowly evolving 
scalar field named quintessence [21]- [IQ], both with a equation of state (EOS) parameter 
w = Pde/pde varying with the redshift. In most of the cosmological models considered 
in the literature, it is generally assumed that dark matter (DM) and DE interact only 
gravitationally. However, it has been argued that, in the absence of an underlying 
symmetry suppressing a possible coupling between dark matter and dark energy, a 
phenomenological interaction between them is not only possible but necessary |43j-[H]. 
In addition to, it has been speculated that an interaction between dark components is 
able to alleviate or inclusive to solve the problem of the Cosmic Coincidence |43j-|48j. 
As a consequence, a large body of work dealing with such a possibility has been explored 
in the literature [l9]-[50]. 

However, the existence or not of some class of interaction between dark components 
is to be discerned observationally. To this respect, constraints on the strength of such 
interaction have been put using different observations [51j-[89j. 

Recently, it has been suggested that an interacting term Q{z) dependent of the 
redshift crosses the noninteracting line Q{z) = [88]- [89]. In [88], this conclusion have 
been obtained using observational data samples in the range z G [0, 1.8] in order to 
fit a scenario in which the whole redshift range is divided into a determined numbers 
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of bins and the interaction function is set to be a constant in each bin. They found 
an oscillatory behavior of the interaction function Q{z) changing its sign several times 
during the evolution of the universe. On the other hand, in |89] is reported a crossing of 
the noninteracting line Q{z) = under the assumption that the interacting term Q{z) is 
a linearly dependent interacting function of the scale factor with two free parameters to 
be estimated. They found a crossing from negative values at the past (energy transfers 
from dark matter to dark energy) to positive values at the present (energy transfers 
from dark energy to dark matter) at2;~0.2 — 0.3. 

In order to shed light on the question of if an interaction term can solved the 
The Cosmic Coincidence problem or if in this scenario such crossing really exists, in a 
previous article jl] we propose the reconstruction of a quite general nongravitational 
interaction Q between dark components which we introduce phenomenologically in 
the equations of motion. This reconstruction of Q as a function of the redshift 
was done in terms of an expansion of Chebyshev Polynomials which constitute a 
complete orthonormal basis on the finite interval [-1,1] and have the nice property 
to be the minimax approximating polynomial (this technique has been applied to 
the reconstruction of the DE potential in |9n]-[91])- To this respect, we did the 
reconstruction of the interaction function using the data from the Union2 la Supernova 
(SNIa) data test (557 data) [3]. In this article, we use a set of data samples covering 
the wide range of redshift z G [0, 1089] in order to reconstruct the interaction function 
Q from late to early times. 

A summary of the paper is as follows. In the second section, we introduce the 
general formalism and the equations of motion representing a DE fluid interacting with 
a DM fluid in Einstein gravity. In the third section, we write the cosmological equations 
of motion for the interacting dark fluids in a Friedmann-Robertson- Walker (FRW) 
spacetime. We continue in the section fourth with the introduction and development 
of the method for reconstructing the interaction function in terms of an expansion 
in the base of Chebyshev polynomials. The fifth section is dedicated to describe the 
parametric bayesian statistical method used in the reconstruction of the interaction 
function together with the chi-square test and the priors used on the free parameters 
of each model. Additionally, we present the observations used in this article in order to 
fit the free parameters of the models considered. These observations consist in the data 
from the Union2 la Supernova (SNIa) data test (557 data) [3j, the BAO from SDSS 
DR7 [8], the CMB from 7- year WMAP experiment ^Q\, the Hubble expansion rate (15 
data) [93]- [91], [95] and the X-ray gas mass fraction (42 data) [96] . 

Then, the sixth section is dedicated to the reconstruction of the interaction function 
and, to put constraints on the parameters of each model and on the coefficients of the 
Chebyshev expansion Then we give a discussion of the results of the reconstruction and 
the best estimated values of the parameters of every model fitting the observations. 
Finally, in the last section we present a brief resume of our principal results and our 
conclusions. 
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2. General equations of motion for dark energy interacting with dark 
matter. 

We assume an universe formed by four components: the baryonic matter fluid (6), 
the radiation fluid (r), the dark matter fluid (DM) and the dark energy fluid {DE). 
Moreover all these constituents are interacting gravitationally and additionally only the 
dark components interact nongravitationally through an energy exchange between them 
mediated by the interaction term deflned below. 

The gravitational equations of motion are the Einstein fleld equations 



G^y = 8ttG [T^^ + T^^ + T^^ + T^^] , (1) 
whereas that the equations of motion for each fluid are 

VXu = 0, (2) 

v-t;, = 0, (3) 

V%^/ = -F,, (4) 

V%^/ = F,, (5) 

where the respective energy-momentum tensor for the fluid i is deflned as {i = 
b, r,DM,DE), 

= Pi ^ii^y + + UtMUu)Pi, (6) 

here is the velocity of the fluids (assumed to be the same for each one) where as pi 
and Pi are respectively the density and pressure of the fluid i measured by an observer 
with velocity u'^. F^ is the cuadrivector of interaction between dark components and 
its form is not known a priori because in general we do not have fundamental theory, 
in case of existing, to predict its structure. We project the equations ([S])-© in a part 
parallel to the velocity u^, 

u^VXu = 0, (7) 

u^VXu = 0, (8) 

u^VXu'' = (9) 

u^VXJ" = ^'F,^ (10) 

and in other part orthogonal to the velocity using the projector hp^ = gp^ + UfjU^ acting 
on the hypersurface orthogonal to the velocity u^, 

h^f'vXu = 0, (11) 

h^^VXu = 0, (12) 
ht^f^^-Tl^M ^ _h^PF^, (13) 
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using ([6]) in (l71)-(IT0l) we obtain the mass energy conservation equations for each fluid, 
M^V^Pfe + (Pfc + Pfe) VX = 0, (15) 

U^V^Pr + {Pr + Pr) = 0, (16) 

u^^V^Pdm + {pDM + Pdm) V^u^ = W'F^, (17) 

u^V^PDE + {PDE + Pde) = -u^F^- (18) 

On the other hand, introducing ([6]) in f|TT|) - f|T^ it permits to have the Euler equations 
for every fluid, 

h^^V^P, + (pb + P,) u^Vy = 0, (19) 

hf^^Vf^Pr + {pr + Pr) u'^V ^ = 0, (20) 

h^^V.PoM + {pDM + Pdm) u^V.u^ = -h^^F^, (21) 

h^^V^PoE + {pDE + Pde) u^V^ = F^. (22) 

Finally we closed the system of equations assuming the following state equations for the 
respectively baryonic, dark matter, radiation components, 

Pb = 0, (23) 
Pdm = 0, (24) 

Pr =\pr, (25) 

while for the dark energy we assume a state equation with constant parameter 

Pde = wpDE- (26) 

3. Cosmological Equations of motion for dark energy interacting with dark 
matter. 

We assumed that the background metric is described by the flat Friedmann-Robertson- 
Walker (FRW) metric written in comoving coordinates as supported by the anisotropies 
of the cosmic microwave background (CMB) radiation measured by the WMAP 
experiment [9] 

ds^ = -dt" + a^{t) {dr^ + r^dn^) , (27) 

where a{t) is the scale factor and t is the cosmic time. In these coordinates we choose 
for the normalized velocity, 

m'^ = (1,0,0,0), (28) 

and therefore we have, 

= 3 ^ = 3i/, (29) 



w 



= 0, (30) 
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where H is the Hubble parameter and the point means derivative respect to the cosmic 
time. In congruence with the symmetries of spatial isotropy and homogeneity of the 
FRW spacetime, the densities and pressures of the fluids are depending only of the 
cosmic time, Pi(t), Pi(t), and at the same time, the parallel and orthogonal components 
of the cuadrivector of interaction with respect to the velocity are respectively, 

M^F^ =Q{a), (31) 
h^^F^ = 0, (32) 

where Q{a) is known as the interaction function depending on the scale factor. The 
introduction of the state equations ( l23l) - (!26ll . the metric (|27ll and the expressions (l28l) - 
f p2|) in the equations of mass energy conservation for the fluids f lT5|) -f lT8|) produces, 

Pb + 3Hpb = 0, (33) 

Pr + AHpr = 0, (34) 

Pdm + SHpDM = Q, (35) 

Pde + 3{1 + w) HpDE = -Q. (36) 

On the other hand, the Euler equations (fT9|) -( l22|) are satisfied identically and do not 
produce any new equation. From the Einstein equation ([T]) we complete the equations 
of motion with the first Friedmann equation. 

Its convenient to define the following dimensionless density parameters Q*, for i = 
b, r, DM, DE, as the energy densities normalized by the critical density at the actual 
epoch, 

- -^, (38) 

Pcrit 

and the corresponding dimensionless density parameters at the present, 

^° = (39) 

Pcrit 

where p^^.^^ = SHq/StiG is the critical density today and Hq is the Hubble constant. 
Solving f l33|) and fl34|) in terms of the redshift z, defined as a = 1/(1 + z), we obtain the 
known solutions for the baryonic matter and radiation density parameters respectively: 

ntiz) = n^ii + zf, (40) 
n:.{z) = n',{i + z)\ (41) 

The energy conservation equations (l35l) and (l36i) for both dark components are rewritten 
in terms of the redshift as: 

dpDM 3 Q{z) 

dz I + z [1 + z) ■ H[z) 

(^Pde 3(1 + w) Q{z) 
dz I + z [1 + z) ■ H[z) 
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Phenomenologically, we choose to describe the interaction between the two dark fluids 
as an exchange of energy at a rate proportional to the Hubble parameter: 



Q{z)^pl,,.{l + zf.Hiz)-Uz). (44) 

The term p'^^-^ ■ (1 + zY has been introduced by convenience in order to mimic a rate 
proportional to the behavior of a matter density without interaction. Let be note that 
the dimensionless interaction function Iq{z) depends of the redshift and it will be the 
function to be reconstructed. With the help of (HI]), we rewrite the equations for the 
dark fluids (US])- (03]) as, 



DM 



dz 



dn 



DE 



Q 

l + z 

3(l + w 



dz 



U, = -a + zf-lQ{z), 



(45) 
(46) 



4. General Reconstruction of the interaction using Chebyshev polynomials. 



We do the parametrization of the dimensionless coupling Iq{z) in terms of the Chebyshev 
polynomials, which form a complete set of orthonormal functions on the interval [—1, 1]. 
They also have the property to be the minimax approximating polynomial, which means 
that has the smallest maximum deviation from the true function at any given order [90] - 
[5T] . Without loss of generality, we can then expand the coupling Iq(;z) in the redshift 
representation as: 

N 

1q(z)^^A„-T„(^), (47) 

n=0 

where Tn{z) denotes the Chebyshev polynomials of order n with n G [0, A^] and 
a positive integer. The coefficients of the polynomial expansion A„ are real free 
dimensionless parameters. Then the interaction function can be rewritten as 



N 



0° 

rent 



[l + zf ■H{z)-Y,^n-Tn{z). 



(48) 



n=0 



We introduce ( H71) in (I45l) - (l46l) and integrate both equations obtaining the solutions. 



DM\ 



[l + z) 



N 



^l-DM 



5^A„ ■Kn{x,0) 



n=0 



^DEiz) 



:i + z 



,3(l+to) 



N 



0° 



W) 



n=0 



where we have defined the integrals 

Tnix) 



Kn{x,w) 



1 (a + 6£)(i+3«') 



dx 



(49) 
(50) 

(51) 



Dark energy interacting with dark matter 



8 



and the quantities, 



x = ^-l, (52) 



^max 

Z. 



a = l + (53) 



b = (54) 

here Zmax is the maximum redshift at which observations are available so that x E [—1,1] 
and \Tn{x)\ < 1, for all n e [0,N]. 

Finally, using the solutions fj40|) - PT|) and (H9|) - (!50|) we rewrite the Friedmann equation 
(EZD as 

(z) = Hi [fi°(l + zf + fi°(l + zf + fiW(^) + ^de{^)\ ■ (55) 

The Hubble parameter depends of the parameters {Hq, VL^, fi^, fi^jj^^, ^^^e^ w) and the 
dimensionless coefficients A^. However one of the parameters depends of the others due 
to the Friedmann equation evaluated at the present, 

n'^E = i-nl~nl-nUj. (56) 

At the end, for the reconstruction, we have the five parameters {Hq, w) 
and the dimensionless coefficients A„. 

To do a general reconstruction in ( I49l) -( l50l) we must take — t- oo and to obtain 
the solutions in a closed form. The details of the calculation of the integrals Kn{x,w) 
in the right hand side of (l49l) -( |50l) are shown in detail in the Appendix A which shows 
the closed forms (A.9)-(A.10) for the integrals with odd and even integer n subindex, 
and valid for w ^ n/3, where n > (see ref. [IJ). Finally, we point out the formulas we 
use for the reconstruction of other important cosmological properties of the universe: 

• The age of the universe at redshift z in |92j : 

• The deceleration parameter: 

5. Current observational data and cosmological constraint. 

To simplify our analysis, we reconstruct the coupling function Iq(-2) to different orders 
(A^ = 1,2,3,4), up to order A^ = 4. The details of this reconstruction are described in 
the Appendix B. (see ref. \T\). In order to do it, we test our cosmological models using 
the observational data currently available. In this section, we then describe how we use 
these data. 
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5.1. Type la supernovae. 

We test and constrain the coupling function Iq(-2) using the "Union2" SNe la data set 
from "The Supernova Cosmology Project" (SCP) composed by 557 type la supernovae 
[3]. As it is usual, we use the definition of luminosity distance (11 (see ref. [2]) in a flat 
cosmology, 

d4..X)=c(l + .)y^ (59) 

where H{z,X.) is the Hubble parameter, i.e., the expression (^S^, "c" is the speed of 
light given in units of km/sec and X represents the parameters of the model, 

X = {Ho, nl qIm, Ai, a^v). (go) 

The theoretical distance moduli for the k-th supernova with redshift Zk is defined as 



/i*^(;zfc, X) = m{z) - M = 5 log 



10 



dL{zk,X.) 
Mpc 



+ 25, (61) 



where m and M are the apparent and absolute magnitudes of the SNe la respectively, 
and the superscript "th" stands for "theoretical" . We construct the statistical Xsn 
function as 

xlN(X).f i^MzM!, (62, 

k=l 

where /i^ is the observational distance moduli for the k-th supernova, al is the variance 
of the measurement. With this function xIn' construct the probability density 
function (pdf) as 

pdfsN(X) = Ai-e->^iNW/2, (63) 
where Ai is a integration constant. 

5.2. Baryon acoustic oscillations 

The baryon acoustic oscillations (BAO) are detected in the clustering of the combined 
2dFGRS and SDSS main galaxy samples, and measure the distance-redshift relation at 
z = 0.2. Additionally, BAO in the clustering of the SDSS luminous red galaxies measure 
the distance-redshift relation aX z = 0.35. The observed scale of the BAO calculated 
from these samples, as well as from the combined sample, are jointly analyzed using 
estimates of the correlated errors to constrain the form of the distance measure D^{z) 

where X represents the parameters of the model, see equation (160|) . The peak position 
of the BAO depends on the ratio of Dy{z) to the sound horizon size at the drag epoch 
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The observational data 



z 


dz 


0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



Table 1. Summary of the BAO data set. 



(where baryons were released from photons) z^, which can be obtained by using a fitting 
formula [1]: 

where Qm = ^dm + 

h = 0.313(nM)-°-''' [1 ± 0.607(1]m)°-'''] , (66) 
62 = 0.238((]m)°-'''. (67) 

We use the distance radio d^ aX z = 0.2 and z = 0.35 |1]-[Z] 

"^'■'^^^ = Dy{o.2,^y "^'-''^^^ = DJiolt^y ^^^^ 

where rs{zd, X) is the comoving sound horizon size at the baryon drag epoch. 

1 dz' 
r,{z,yi) = —= " (69) 

Using the data of BAO of the table [1] and the following inverse covariance matrix of 
BAO in [7] 

1 _ ( ±30124 -17227 \ 
^BAO- _i7227 ±86977 ) 

thus, the function of the BAO data is constructed as: 

Xbao(X) = (rff (X) - df^y [CAo],, (rf*'(X), - rff ) , (71) 

where [d^^ — d°^'^) is a column vector formed from the values of theory minus the 
corresponding observational data, with 

jth(^\ robs _ ( c?o.2o(X) - 0.1905 \ 
^^^""^"^^ -|^rfo.35(X)- 0.1097 ) 

and "t" denotes its transpose. With this Xbao function, we construct the probability 
density function (pdf) as 

pdfB^o(X) = A2-e->^BAo(x)/2 (73) 

where A2 is a integration constant. 



Dark energy interacting with dark matter 



11 





302.40 




1.7239 




1090.89 



Table 2. The observational CMB data. 



5.3. Cosmic microwave background (CMB) 
The CMB shift parameter R is provided by 



R(z„X) = Y^fiMH2(l + z,)(DA(z„X)/c), (74) 

where X represents the parameters of the model (see equation fl60|) and the second 
distance ratio is given by: 

1 r^* dz' 

Da{z.,^) = J-—- JTTT^y (75) 

(1 + z^)J H{z',X.) 

The redshift z^: (the decouphng epoch of photons) is obtained using the fitting function 



z, = 1048 [1 + 0.00124(fi°)-°-^38j ^ g.inMY'] , (76) 

where flu = + functions gi and g2 are 

^ 0.0783(l],o)-°-^3« 

1 + 39.5(^0)°-^''' 
0.560 

- l + 21.1(ao)i-«i- ^^^^ 
In addition, the acoustic scale /a describes the ratio Da{z^,) /rg^z^,), and is defined as 

lAiX)^il + z,) " (79) 

where a factor of (1 + z^) arises because Da{Z:^,X) is the proper angular diameter 
distance, whereas rs{z^,, X) is the comoving sound horizon at z^,. The fitting formula of 
rs{z,X.) is given by the equation (!69l) . The seven year WMAP observations [10] give 
the maximum likelihood values according to table El Following [in], the function for 
the CMB data is 

x'cmb(X) = (xf (X) - )* (CcU)., - ' (80) 

where [x*^ — x°^'^) is a column vector formed from the values of theory minus the 
corresponding observational data, with 

/ ^(2,)- 302.40 \ 
xf{X) - xf' = R{z,) - 1.7239 

^z^- 1090.89 y 
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z 


0.0 


0.1 


0.17 


0.27 


0.4 


0.48 


0.88 


0.90 


1.30 


1.43 


1.53 


1.73 


H{z){Kms-^Mpc~^) 


74.2 


69 


83 


77 


95 


97 


90 


117 


168 


177 


140 


202 


la 


±3.6 


±12 


±8 


±14 


±17 


±60 


±40 


±23 


±17 


±18 


±14 


±40 



Table 3. The observational H{z) data. 



"t" denotes its transpose and {C(-^-Q)ij is the inverse covariance matrix. In [TO], the 
inverse covariance matrix is also given as follows 

/ ±2.3050 ±29.6980 -1.3330 \ 
CcMB = ±29.698 ±6825.27 -113.18 . (82) 
y -1.3330 -113.180 ±3.4140 y 

With this function XcmB) construct the probability density function (pdf) as 

pdfcMB(X) = A3-e-^cMB(x)/2 (83) 

where A3 is a integration constant. 
5.4- Hubble expansion rate (H) 

The Hubble parameter can be written as the following form: 

, , 1 (aj Z , , 

= * ■ <«^' 

So,through measuring dz/dt, we can obtain H{z). In [99], jlOO] and [90], the authors 
indicated that it is possible to use absolute ages of passively evolving galaxies to compute 
values of dz/dt. The galactic spectral data used by [90] come from the Gemini Deep 
Survey [101] and the archival data [102j - [T07] . Detailed calculations of dz/dt can be 
found in [HO]- Simon obtained H{z) in the range of [0,1.8]. The twelve observational 
Hubble data from [93], [91] are list in the tableO In addition, in [95], the authors took the 
BAO scale as a standard ruler in the radial direction, obtain three more additional data: 
H{z = 0.24) = 79.69±2.32, H{z = 0.34) = 83.80±2.96 and H{z = 0.43) = 86.45±3.27 
(in units of Km s~^Mpc~^). 

We can use these data to constraint different kinds of dark energy models and 
determine the best fit values of the model parameters by minimizing [109] 

X^(X).f ["■^'"■-)-"°'-'-)]' . (85) 



rth 



where X represents the parameters of the model (see equation (l60l) . H is the theoretical 
value for the Hubble parameter, H°'"^ is the observed value, (j{zi) is the standard 
deviation of the measurement, and the summation is over the 15 observational Hubble 
data points. This test has already been used to constrain several cosmological models 
|110] - |123] . With this function xh. '^^ construct the probability density function (pdf) 
as 

pdfH(X) = A4-e-'^H(x)/2 (gg) 
where A4 is a integration constant. 
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5.5. X-ray gas mass fraction (X-ray) 

According to the X-ray cluster gas mass fraction observations, tlie baryon mass fraction 
in clusters of galaxies can be used to constrain cosmological parameters. The X-ray 
cluster gas mass fraction, fgas, is defined as the radio of the X-ray gas mass to the total 
mass of a cluster, which is a constant and independent of the redshift. In the framework 
of the ACDM reference cosmology, the X-ray cluster gas mass fraction is presented as 

ESI 



where X represents the parameters of the model. The factor K is used to describe 
the combined effects of the residual uncertainties, such as the instrumental calibration 
and certain X-ray modeling issues, and a Gaussian prior for the calibration factor is 
considered byi^' = 1.0±0.1 [96]. The parameter 7 denotes permissible departures from 
the assumption of hydrostatic equilibrium, due to non-thermal pressure support; the 
bias factor b{z) = &o(l + o^bZ) accounts for uncertainties in the cluster depletion factor; 
s{z) = so(l + CisZ) accounts for uncertainties of the baryon mass fraction in stars and 
a Gaussian prior for sq is employed, with sq = 0.16 ± 0.05 [96j and A is the angular 
correction factor, which is caused by the change in angle for the current test model 62500 
in comparison with that of the reference cosmology 025^^*^- 

, _ f^Mil\ „ ~ / H{z,X)DAz,X) \ ^ 

(89) 

here, the index rj is the slope of the fgasi^ /r25oo) data within the radius r25oo, with 
the best-fit average value rj = 0.214 ± 0.022 [96j. And the proper angular diameter 
distance Da{z) is given by the equation fl75l) . Following the method in [96] and [124] 
and adopting the updated 42 observational fgas data in [96], the best fit values of the 
model parameters for the X-ray cluster gas mass fraction analysis are determined by 
minimizing. 

With this function Xxray' construct the probability density function (pdf) as 

pdfx^^^(X)=A5-e-^W(^)/2 (91) 

where A5 is a integration constant. We have considered the following values to do the 
fit 

K = +1.0, 7 = +0.9980, So = +0.16, 
as = +0.150, 60 = +0.787, at, = -0.080. 
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Therefore, using the equations (15^ . ( 175]) . and (^1]) we build the total 

probability density function as 

Pdf = pdfsN ■ pdfBAO • pdfcMB ■ Pdfn " pdfx^ay' (92) 
We can then rewrite the above equation as: 

Pdf(X) = A■e-^'W/^ (93) 
from here we define the total function Xt ^s, 

= X^SN + X^BAO + X^CMB + X^H + X^Xray • (94) 

Thus the constraints on cosmological models from a combination of the above discussed 
observational datasets can be obtained by minimizing the equation ( IMI) . 

5.6. Priors on the total probability density function fPdf). 

In the models I, II and III shown in the Table HI we marginalize the parameters 

Y = (f^o, ^i)M) ^fei ^r) the pdf (1921) choosing priors on them. In order to it, we 
must compute the following integration, 

/•oo /"OO /"OO /"OO 

Pdf(v)=/ / / / Pdi{'K)Pdi{Y)dHodn'jjMdnldn'^,, (95) 

Jo Jo Jo Jo 

where V = {w, Ai, X^) represents the nonmarginalized parameters, Pdf (X) is given 
by (|92l) and Pdf(Y) is the prior probability distribution function for the parameters 
{Hq, fi2)M5 ^r) which are chosen as Dirac delta priors around the specific values 

Y = {Hq, ^%My ^r) measured by some other independent observations, 

Pdf (Y) = s{Ho - Ho) ■ 5(^^?)A/ - ^Ui) ■ - r2°) ■ 5{n^. - n^.). (96) 

Introducing (196!) in ( (!95l) it produces, 

Pdf(V) = A-e-*'(^)/^ (97) 

where we have defined a new function depending only on the parameters V = 
{w,Xi,...,Xn) as, 

X\y) ^ E ^ ^. (98) 



The specific values chosen for the Dirac delta priors are, 

• Hq = 70.4 (km/s)Mpc~^ as suggested by the observations in [10]. 
. = 0.227, 

• nl = 0.0456, 

• n° = 8.42095 X 10"^ 
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Models 



Models 


i 'DM 


EOS parameter w 


Interaction function 


Model 1 
Model 11 
Model 111 
Model IV 


0.227 (fixed) 
0.227 (fixed) 
0.227 (fixed) 
free parameter 


constant 
-1 

constant 
constant 


Iq(^)=0 
Iq(^) 7^ 
Iq(^) 7^ 



Table 4. Summary of the models studied in this work. In the models II, III and IV 
the interaction function Iq(z) will be reconstructed. Additionally, in the model IV the 
parameter ri^JM estimated. 



Once constructed the function f l94|) . we numerically minimize it to compute 
the ''best estimates" for the free parameters of the model: V = {w, Xi, X^). The 
minimum value of the function gives the best estimated values of V and measures 
the goodness-of-fit of the model to data. For the Model IV, we leave too the parameter 
free to vary and estimated it from the minimization of the function. In this case, 
the parameters to be marginalized are Y = {Hq, Q^, Q^). Then, the marginalization will 
be as, 

/•oo /"OO /»oo 

Pdf (V) = / / / Pdi{x.)Pdi{Y) dHo dn'^^dn'^^ (99) 

Jo Jo Jo 

where now V = {w,Q^ijj^.j, Xi, Xn) represents the nonmarginalized parameters to be 
estimated, Pdf(X) is given by fl92|) and Pdf(Y) is the prior probability distribution 
function for the parameters {Ho, Q^, f2°) which are chosen as Dirac delta priors around 
the specific values Y = (Hq, fl^, fl^) given above. In the models 11, 111 and IV the 
interaction function 1q(-2) will be reconstructed up to order = 4 in the expansion in 
terms of Chebyshev polynomials. 



6. Results of the reconstruction of the interaction function. 

Now, we present the results of the fit of the models listed in the Table S] with the 
"Union2" SNe la data set [3], the baryon acoustic oscillation (BAO), the cosmic 
microwave background (CMB) data from 7-year WMAP, the observational Hubble data 
and the cluster X-ray gas mass fraction, using the priors described in the Section 15.61 
For the noninteracting model 1, the only free parameter to be estimated is 6 = {w}, 
whilst for the interacting models II, 111 and IV the free parameters are 9 = {Aq, Ajv}, 
9 = {w, Xo, Xn} and 9 = {w,Q%j^,j, Xq, Xn} respectively, where A^ is taking the 
values A^ = 1,2,3,4. In every case, we obtain the best fitted parameters and the 
corresponding xLn- 
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Model I: w = constant. 

Best estimate for the EOS parameter w. 



Errors 


±1(7 


±2(T 


u 


i.UiOU_o.0163 


i.UiOU_o.0334 



Table 5. The best estimates of the dark energy EOS parameter w for the Model I. It 
was computed through a Bayesian statistical analysis using SNela + BAO + CMB + 
H + X-ray data sets giving xlim = 719.2342 



Age of the universe 



Model I 


±1(T (Gyr) 


±2(T (Gyr) 


Age 


io.ou /y_Q 0201 


1 S 8fi7q+° °3'^^ 

iO.OU ( J_o.0400 



Table 6. Best estimates for the age of the universe and their errors at la and 2a for 
the model I using SNela + BAO + CMB + H + X-ray data sets. 



Model II: w = -1. 



Best estimates for the parameters A, 



An 


N = 


N = 1 


A = 2 


A = 3 


A = 4 


Ao 


-1.00 X 10-5 


-1.0296 X 10-5 


-1.0108 X 10-5 


-1.0105 X 10-5 


-1.0120 X 10-5 


Ai 


0.0 


2.0027 X 10-^ 


2.0024 X 10-^ 


2.0014 X 10-^ 


2.0006 X 10-^ 


A2 


0.0 


0.0 


1.9990 X 10-6 


1.9935 X 10-6 


1.9918 X 10-6 


As 


0.0 


0.0 


0.0 


1.00 X 10-" 


1.018 X 10-1^ 


A4 


0.0 


0.0 


0.0 


0.0 


1.00 X 10-16 


■ 

A. mm 


719.8193 


719.8154 


719.8134 


719.8131 


719.8127 



Table 7. Summary of the best estimates of the dimensionless coefficients A„ of the 
polynomial expansion of Iq(z) for the Model II, corresponding to a interacting dark 
energy EOS parameter w = —1. They were computed through a Bayesian statistical 
analysis using the SNela-t-BAO+CMB+H+X-ray data sets. The number N in the top 
of every column indicates the maximum number of Chebyshev polynomials used in the 
expansion of Iq{z) starting from iV = 1 to iV = 4. From the Figure |3] to Figure E] show 
the reconstruction of several cosmological variables using these best fitted values. 



Model II: w = -1. 



Errors 


±la 


±2(7 


Ao 


-i.0108 X iO _o.oio4xio-5 


1 ninS V 10-5+00282x10-^ 
-i.0108 X iO _o.0244xlO-5 


Ai 


1 nno/i 1 ri-4+0-0253xl0-* 
+2.0024 X 10 '_o.0223xlO-4 


1 nno^ -,^.-4+00408x10-* 

+2.0024 X 10 ''.o.osssxio-* 


A2 


, 1 nnnn t n-6+0-0011 X IQ-" 

+1.9990 X 10 6_oooi7^io_6 


, 1 nnnn 1^-6+00028x10-'' 

+1.9990 X 10 6_ooo47^io_6 



Table 8. Summary of the la and 2a errors of the best estimates for N = 2. 
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Age of the universe 



Model II 


±la (Gyr) 


±2a (Gyr) 


Age 


iO.OdUZO_Q QQQQ25 


±O.OOUZO_Q QQQQ39 



Table 9. Best estimates for the age of the universe and their errors at la and 2a for 
the model II using SNela + BAO + CMB + 11 + X-ray data sets. 
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0.008 
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0.006 










0.00015 








0.004 










0.0001 








0.002 










5e-05 




























0.002 










-5e-05 








0.004 










-0.0001 








0.006 






la 

2a 




-0.00015 






■o 


0.008 




Union2 SNIa+BAO+CMB+H+X-ray 


Minimum + 




-0.0002 




Unlon2 SNIa-tBAO+CMB+H+X-ray 


2o 

Minimum + 



-0.0015 -0.001 -0.0005 0.0005 0.001 0.0016 -0.008 -0.006 -0.004 -0.002 0.002 0.004 0.006 0.008 



Figure 1. Contours correspond to la, 2a confidence levels constrained from joint 
SNela + BAO + CMB + 11 + X-ray data analysis for the marginalized probability 
densities of the model II, using an expansion in terms of the first N = 2 Chebyshev 
polynomials in the (Aq - Ai) and (Ai - A2) planes. It is clear that before marginalization, 
we have three free parameters (Aq, Ai, A2). In every figure, we marginalized on one of 
the parameters. 



Model III: w = constant. 
Best estimates for the parameters and w. 



An 


N = 


n = 1 


n = 2 


n = 3 


n = 4 


Ao 


-1.00 X 10"^ 


-1.0040 X 10"^ 


-1.0015 X 10^^ 


-1.0001 X 10-5 


-1.00 X 10-5 


Ai 


0.0 


2.0021 X 10-^ 


2.0040 X 10-^ 


2.0120 X 10-^ 


2.0390 X 10-^ 


A2 


0.0 


0.0 


9.9980 X 10-^ 


9.9970 X 10-^ 


9.9962 X 10"^ 


A3 


0.0 


0.0 


0.0 


1.0010 X 10-^ 


1.0020 X 10-1° 


A4 


0.0 


0.0 


0.0 


0.0 


1.00 X 10-12 




-1.0089 


-1.0088 


-1.0086 


-1.0089 


-1.0092 


■ 

Am«n 


719.2984 


719.2861 


719.1833 


719.1779 


719.1731 



Table 10. Summary of the best estimates of the dimensionless coefficients An of the 
polynomial expansion of Iq(2;) for the Model III, corresponding to a interacting dark 
energy EOS parameter w = constant and fi^jM = 0.227. They were computed through 
a Bayesian statistical analysis using the SNela+BAO+CMB+H+X-ray data sets. The 
number N in the top of every column indicates the maximum number of Chebyshev 
polynomials used in the expansion of Iq(z) starting from = 1 to A'^ = 4. From the 
Figure H] to Figure [8] show the reconstruction of several cosmological variables using 
these best fitted values. 
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Model III: w =constant. 



Errors 


±la 


±2a 


Ao 
Ai 

A2 

u 


i.UUiO X iU -0.0473x10-5 
+ 2 0040 X 10-4+0-5885x10-4 
-I-Z.UU4U X iU -0.1020x10-4 

+9.9980 X io-^s:Sx;;=: 

-1.0086^°:°°- 


1 001^ X in-5+O-ioooxio-'' 

-i.UUlO X iU -0.1339x10-5 

+2.0040 X io-:s:s^;;- 
+9.9980 X io-»«:i?2j;;:: 

-1.0086«:SSt5 



Table 11. Summary of the la and 2cr errors of the best estimates for N = 2. 
Age of the universe 



Model 111 


±la (Gyr) 


±2a (Gyr) 


Age 


15.86064+°:0337 


15.860641°:°^^° 
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Table 12. Best estimates for the age of the universe and their errors at Icr and 2a for 
the model III using SNela + BAO + CMB + H + X-ray data sets. 



-1.06 
-1.08 



-1.1 
-2e- 



Union2 SNIa-fBAO-fCMB-tH-fX-ray 




Modellll 






Iff 




2a 




Minimum + 





I0 




2o 




Ivlinimum + 






Modei iii 




Union2 SNIa-tBAO-tCMB-tH-fX-ray 





-09 -1.5e-09 -1e-09 -5e-10 



-0.009 -0.008 -0.007 -0.006 -0.005 -0.004 -0.003 -0.002 -0.001 0.001 

'-1 



Figure 2. Contours correspond to la, 2a confidence levels constrained from joint 
SNela + BAO + CMB + H + X-ray data analysis for the marginalized probability 
densities of the model III, using an expansion in terms of the first N — 2 Chebyshev 
polynomials in the (Xq-uj) and (A1-A2) planes. It is clear that before marginalization, 
we have four free parameters (Aq, Ai, A2, w). In every figure, we marginalized on the 
last two remaining parameters. Note that the preferred region for the EOS parameter 
w is the phantom region and LCDM model is in the la error region. 



The Figure H] shows the reconstruction of the dimensionless interaction function 
1q{z) as a function of the redshift for the models 11 (corresponding to a dark energy 
EOS parameter w = —1), 111 (corresponding to a dark energy EOS parameter w = 
constant and fixed) and IV (corresponding to a dark energy EOS parameter w = 
constant and fi%M ^ fr^e parameter to be estimated) respectively. 

The Figure [5] shows the reconstruction of the dark matter and dark energy density 
parameters ^2^4^(2;), ^oe{ 

clS cl function of the redshift for the models 11, 111 and IV 

described above. 
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Model IV: w =constant and =constant 

Best estimates for the parameters A„, w and 





n — 
It — u 


'il 1 

It 1 


n — 9 


It 


It — ^ 


An 


—1 nn X in^^ 


—1 00 X 1 0^^ 


—1 001 X 1 0^^ 


—1 0006 X 10^^ 

X.WWWVJ /N xw 


—1 0009 X 10-^ 


Ai 


0.0 


2.0314 X 10-^ 


2.0331 X 10-* 


2.0385 X 10"^ 


2.0220 X 10-^ 


A2 


0.0 


0.0 


9.9990 X 10-6 


9.9925 X 10-6 


9.9915 X 10-6 


A3 


0.0 


0.0 


0.0 


1.0068 X 10-^ 


1.0125 X 10-1° 


A4 


0.0 


0.0 


0.0 


0.0 


1.00 X 10-12 


w 


-1.0089 


-1.0089 


-1.0088 


-1.0088 


-1.0089 


i 'DM 


+0.2278 


+0.2278 


+0.2278 


+0.2279 


+0.2279 


• 

Am«n 


718.6851 


718.6602 


718.6236 


718.5704 


718.5617 



Table 13. Summary of the best estimates of the dimensionless coefficients A„ of 
the polynomial expansion of Iq{z) for the Model IV, corresponding to a interacting 
dark energy EOS parameter w — constant and ^ parameter to be 

estimated. They were computed through a Bayesian statistical analysis using the 
SNela+BAO+CMB+H+X-ray data sets. The number N in the top of every column 
indicates the maximum number of Chebyshev polynomials used in the expansion of 
Iq(z) starting from A'' = 1 to TV = 4. From the Figure H] to Figure [8] show the 
reconstruction of several cosmological variables using these best fitted values. 



Model IV: w =constant and ^^nnr =constant. 



Errors 


±1(T 


±2(7 


Ao 
Ai 
A2 

i 'DA/ 


1 0010 X 10-5+0-0616x10-'' 
-l.UUiU X lU -0.01015x10-5 

+ 2 331 X 10 -4+0.0756x10-4 
X iU -0.0173x10-4 

+9 9990 X lo-e+o-osgixio-l; 

-t-J.yjJU A ±U -0.0489x10-6 

1 nnss+3-9588xio-* 
-J--U088-0.0658 
1 n OO7S+8-0936X10-5 
+U-^-^'O-0.0056 


-1 00 1 X 10 -5 +0.1122x10-;; 

l.UUiU X iU -0.0434x10-5 
+ 2 0331 X 10-4+0-1887x10-4 
-t-Z.UOOi X lU -0.0681x10-4 

+9 9 9 90 X 10 -6+0-1550x10-6 

-hy.JJJU X lU -0.1073x10-6 

1 nnfis+i-8i09xio-^ 

in 0070+2. 58411x10-4 
-t-U.ZZ/ O-0.0079 



Table 14. Summary of the la and 2a errors of the best estimates for N ~ 2. 



Age of the universe 



Model IV 


±la (Gyr) 


±2(7 (Gyr) 


Age 


1 5 851 5+°-!^°^ 

lO.OOlO-o.ooiO 


1 5 851 5+°-i^°^ 

lO.OOiO-0.0032 



Table 15. Best estimates for the age of the universe and their errors at Icr and 2a for 
the model IV using SNela + BAO + CMB + 11 + X-ray data sets. 



The Figure |6] shows the reconstruction of the deceleration parameter q{z) as a 
function of the redshift for the models 11, 111 and IV respectively. 

The Figure [7] shows the reconstruction of the age of the universe HQt{ clS db function 
of the redshift for the models 11, 111 and IV described respectively. 

The Figure E] shows the reconstruction of the rate between dark density parameters 
^d£;(-2)/^dm(^) ^ function of the redshift for the models 11, 111 and IV described 
above. 
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la 


0.98 










Minimum + 


-1 






1.02 






1.04 






1.06 


Model IV ' 




1.08 


Unlon2 SNIa+BAO+CMB+H+X-ray 





0.222 0.224 0.226 0.228 0.23 0.232 0.234 0.236 0.238 0.24 0.242 



0.225 
-0.003 



Unlon2 SNIa+BAO+CMB+H+X-ray 






Model IV 




; 












1a 






2a 






Minimum + 



Model IV 


1o 




2a 




Minimum + 


Union2 SNIa+BAO+CMB+H+X-ray 





-0.0003 

-0.007 -0.006 -0.005 -0.004 -0.003 -0.002 -0.001 



Figure 3. Contours correspond to la, 2a confidence levels constrained from joint 
SNela + BAO + CMB + H + X-ray data analysis for the marginalized probability 
densities of the model IV, using an expansion in terms of the first N = 2 Chebyshev 
polynomials in the {floM - w), (Ao - ^dm) and (Ai - A2) planes. It is clear that before 
marginalization, we have five free parameters (Aq, Ai, A2, uj, il/^jv/)- In every figure, we 
marginalized on the last three remaining parameters. Note that the preferred region 
for the EOS parameter w is the phantom region. Our result is consistent with the 
LCDM model in the la error region. 



All these Figures show the superposition of the best estimates for every cosmological 
variable in terms of the parameters A„ corresponding to the coefficients of the polynomial 
expansion (H7|) ranging from iV = 1 to 4. The Tables [TJ [10] and [13], show the best 
fitted parameters and the minimum of the function Xmin models II, III and IV 

respectively. From these tables, we can note the fast convergence of the best estimates 
when the numbers of parameters is increased in the expansion P7|) . 



The Figure [9] shows the reconstruction of the dimensionless interaction function 
Iq(2;) as a function of the redshift, using the best estimates for N = 2 and their 
confidence intervals at la and 2a for the models II (corresponding to a dark energy 
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0.00036 
0.0003 
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0.0002 
0.00015 
0.0001 
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-5e-05 

0.0004 
0.00035 
0.0003 
0.00025 
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0.00015 
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6e-06 


-6e-06 



Model 11 




Union2 SNIa+BAO+CMB+H+X-ray 


Best fit 




h, 












X^,A.^,A.2-^3 




Xq,X^,X2'^3'^4 




D 0.2 0.4 0.6 0.8 1 


1.2 1.4 1.6 1 


z 




Model IV 




Unlon2SNIa+BAO+CMB+H+X-ray 


Best fit 
















^0'^l'^2'^3 




XQ,'k-^,X2'^3'K 




D 0.2 0.4 0.6 0.8 1 


1.2 1.4 1.6 1 



0.0004 
0.00035 
0.0003 
0.00025 ; 
0.0002 - 
0.00015 
0.0001 
5e-05 


-5e-05 



Union2 SNIa+BAO+CMB+H+X-ray 



Figure 4. Reconstruction for tlie dimensionless interaction function Iq(z) as a 
function of the redshift for the models II (left above panel), III (right above panel) and 
IV (left below panel) corresponding to a dark energy equation of state parameter (II) 
w = —1, (III) w = constant (both with ri^JM — 0.227), and (IV) w — constant, ri^JM — 
constant, respectively. The curves with different colors show the best estimates using 
the expansion of Iq(z) in terms of the parameters A„ corresponding to the Chebyshev 
polynomial expansion ranging from = 1 to 4. Note the fast convergence of the 
curves when the number of polynomials N involved in the expansion increases. The 
reconstruction is derived from the best estimation obtained from the combination of 
SNela + BAO + CMB + H + X-ray data sets. Note that the best estimated values of 
the strength of the interaction cross marginally the noninteracting line Iq(^) = only 
at the present changing sign from positive values at the past (energy transfers from 
dark energy to dark matter) to negative values almost at the present (energy transfers 
from dark matter to dark energy). 



EOS parameter w = —1), III (corresponding to a dark energy EOS parameter w = 
constant and fixed) and IV (corresponding to a dark energy EOS parameter w = 
constant and as a free parameter to be estimated) respectively. 

The Figure do] shows the reconstruction of the dark matter and dark energy density 
parameters fl*^j^j{z), Q*jj^{z) as a function of the redshift, using the best estimates for 
N = 2 and their confidence intervals at Icr and 2a for the models II, III and IV described 
above. 
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Figure 5. Same explanation as Figure [9] but now for the reconstruction of the dark 
matter and dark energy density parameters, r2^j^^(z), r2^^(z), as a function of the 
redshift for the model II (left above panel), III (right above panel) and IV (left below 
panel) respectively. Note that the density parameter of dark energy is definite positive 
for all the range of redshift considered in the reconstruction. 



The Figure [TT] which shows the reconstruction and the behavior of the deceleration 
parameter q{^z) as a function of the redshift and their confidence intervals at Icr and la 
for the models II, III and IV respectively. 

The Figure [12] shows the reconstruction of the age of the universe HQt{z) as a 
function of the redshift and their confidence intervals at la and 2a for the models II, 
III and IV already described respectively. 

The Figure [T3] shows the reconstruction of the rate between dark density parameters 
^d£;(^)/^dm(^) ^ ^ function of the redshift and their confidence intervals at la and 
2a for the models II, III and IV already mentioned. 

These last Figures show the superposition of the best estimates for every 
cosmological variable and their confidence intervals at la and 2a. The Tables [HI [ID 
and [TU show the best fitted parameters with their respective la and 2a. 
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Figure 6. Similar explanation as Figure |4] for the reconstruction of the deceleration 
parameter q{z) as a function of the redshift for the model II (left above panel), III 
(right above panel) and IV (left below panel) respectively. At the same time, they are 
compared with the corresponding curve for the LCDM (Lambda Cold Dark Matter) 
model. 



From Figures m to [9] we notice that, for all interacting models, the best estimates for 
the interaction function Iq(-2) cross marginally the noninteracting line Iq(-z) = during 
the present cosmological evolution (at around z ^ 0.09) changing sign from positive 
values at the past (energy transfers from dark energy to dark matter) to negative values 
at the present (energy transfers from dark matter to dark energy). However, taking in 
account the errors corresponding to the fit using three parameters (A^ = 2), we see that 
within the la and 2cr errors, it exists the possibility of crossing of the noninteracting line 
in the recent past at around the range z G (0.08,0.12). Crossings of the noninteracting 
line Q{z) =0 have been recently reported at the references [HE] (with an interacting 
term Q{z) proportional to the Hubble parameter) and jl]. The direction of the change 
found in our work is in the same direction to the results published by these references 
where a crossing (from positive values at the past to negative values at the present was 
found at 2 ~ — 0.12. On the other hand, we did not find the oscillatory behavior of the 
interaction function found by Cai and Su [88] who, using observational data samples in 
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Figure 7. Same as Figure 21 shows the reconstruction of the age of the universe 
Hot{z) as a function of the redshift for the model II (left above panel), III (right above 
panel) and IV (left below panel) respectively. Here, They also compared with the 
corresponding curve for the LCDM (Lambda Cold Dark Matter) model. 



the range z G [0, 1.8], fitted a scheme in which the whole redshift range is divided into 
a determined numbers of bins and the interaction function set to be a constant in each 
bin. 

From Figures [5] to Figures ITOl show that, for all the interacting models studied in 
this work, the best estimates for the dark energy density parameter Q*jj^{z) become to 
be definite positive at all the range of redshifts considered in the data samples. However, 
this statement is conclusive because within the Icr and 2cr errors for the fit with three 
parameters (A^ = 2), the fl}j^{z) becomes to be positive in all the range of redshifts 
considered (remember that we are using five data sets and then the space parameters 
permitted is remarkably reduced). According to the Model 4, note also that the behavior 
of dark matter fl*^j^{z) density parameter, for la and 2a errors exist the possibility of 
being smaller compared to the noninteracting model (LCDM model). 

From Figures [6] to Figure [HI show that, for all models, a transition from a 
deceleration era at early times dominated by the dark and baryonic matter density to an 
acceleration era at late times corresponding to the present domain of the dark energy 
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Figure 8. Superposition of the best estimated curves for the rate between dark density 
parameters r2^^(2;)/r2^^,j(z) for the model II (left above panel), III (right above panel), 
IV (left below panel). In the above figures, the different colored curves show the best 
estimates using the expansion in terms of the first N = 1, 2, 3, 4 Chebyshev polynomials 
respectively. They also compared with the corresponding curve for the LCDM (Lambda 
Cold Dark Matter) model. 



density. This transition deceleration-acceleration to take place at redshift at around 
z ~ 0.7 (it was also obtained by |129j ). while that the rhythm of current universe at 
around go ~ 0.6. It is clear for our constraint, that the universe tends to an early time 
to accelerate and a milder expansion rhythm at present. 

We can see from the figure [5] to [10] that, at the present, the dark energy density 
parameter becomes Q%^{z) ^0.7, which is sufficiently large to generate a non negligible 
dimensionless interacting term of the order of Iq ~ — 10~^, as is shown in the Figures H] 
andO In fact, in these same figures we can appreciate that in the interval of redshifts 
z G [0,1090.89], the dimensionless interaction is in the range Iq G [— 10~^,50] for the 
Model II, Iq G [-10-^35] for the Model III and Iq G [-10-^25] for the Model IV 
respectively, corresponding at 2a error. The order of magnitude of this interaction is 
in agreement with the local constraints put on the strength of a constant dimensionless 
interaction derived from the fit to a data sample of virial masses of relaxed galaxies 
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Figure 9. Comparison of tlie best fitted reconstructed dimensionless interaction 
function Iq(z) and its errors, using the expansion of Iq(^) in terms of the parameters 
Ao, Ai, A2 and the corresponding Chebyshev polynomials for the models II (left above 
panel), III (right above panel) and IV (left below panel) corresponding to a dark energy 
equation of state parameter (II) w — (III) w = constant (with f^£)A/ = 0.227) and 
(IV) w = constant and ri^JA/ ~ constant, respectively. Red lines show best fitted 
reconstructed results, while green lines and blue lines show reconstructed errors within 
the Icr and 21a confidence level errors, obtained from a combination of SNela + BAO 
+ CMB + H + X-ray dataset. Otherwise, note that within the la and 2a errors it 
could be the possibility that the crossing of the noninteracting Iq(2;) = line happens 
before the present. 
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Figure 10. Same as Figure El comparison of tlie best fitted reconstruction and its 
errors of the dark matter r2^^(z) and dark energy iljj^(z) density parameters, in 
terms of the parameters Aq, Ai, A2 and the corresponding Chebyshev polynomials 
for the models II (above panel), III (right above panel) and IV (left below 
panel), respectively. Red lines show best fitted reconstructed results, while green lines 
and blue lines show reconstructed errors within the Icr and 2\a confidence level errors, 
determined from a combination of SNela + BAO + CMB + H + X-ray dataset. In 
addition, note that within the \a and 2cr errors, our best fitted results and its errors 
are consistent with the constraints on f2Jj^^(z) and f2J)^(z) in the LCDM model. 
Furthermore, from the last Figure (right below panel) at early time, the Itr and 2<t 
constraints for ^^^^(z) are above of the corresponding for the LCDM model, this 
considerably alleviates the coincidence problem albeit it does not solve it in full. 
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Figure 11. Same as Figure [9] Comparison and evolution of the best fitted 
reconstructed and its errors of the deceleration parameter q{z), in terms of the 
parameters Aq, Ai, A2 and the corresponding Chebyshev polynomials for the models 
II (left above panel), III (right above panel) and IV (left below panel), respectively. 
Red lines show best fitted reconstructed results, while green lines and blue lines show 
reconstructed errors within the la and 21a confidence level errors, determined from a 
combination of SNela + BAO + CMB + H + X-ray dataset. In addition, note that 
within the la and 2a errors, our result and its errors are consistent with the constraints 
on q{z) in the LCDM model. All Figures show a strong evidence of acceleration in the 
recent past which is consistent with a previous study in |126j . Also the strongest 
evidence of acceleration again happens around the redshift z 0.2. The transition 
redshift when the universe underwent the transition from deceleration to acceleration 
is found to be z 0.7 at the la level. Furthermore, our best fit results and its errors at 
la and 2a are in LCDM model. 



clusters obtained using weak lensing, x-ray and optical data |127] . 

On the contrary, a recent study fitting CMB anisotropy data from the seven-year 
Wilkinson Microwave Anisotropy Probe (WMAP) [10] , the BAO distance measurements 
[7\, the Constitution sample of Snia [2] and constraints on the present-day Hubble 
constant, put stronger constraints on the magnitude of such dimensionless strength of 
the order of ^ ^ lO^^ - IQ-^ [T28] . 
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Figure 12. Same as Figure HI Evolution of the best fit reconstructed and its 
errors of the age of the universe Hot{z), in terms of the parameters Aq, Ai, A2 and 
the corresponding Chebyshev polynomials for the models II (left above panel), III 
(right above panel) and IV (left below panel) , respectively. Red lines show best fit 
reconstructed results, while green lines and blue lines show reconstructed errors within 
the la and 2la confidence level errors, determined from a combination of SNela + 
BAO + CMB + H + X-ray dataset. In addition, note that within the la and 2a 
errors, our result and its errors are consistent with the constraints on Hf)t{z) in the 
LCDM model. Furthermore, our best fit results and its errors at la and 2a are in 
LCDM model. 



Another important aspect concerns to study the coincidence problem, then from 
the Figures 151 and fTOl we note that the Models II and III, do not have the possibility of 
alleviating this problem. Their behaviors are in a degree similar to the LCDM model. 
On the contrary, our Model IV, has the possibility of alleviating it at 2o" error. In 
the behavior of dark matter ^^^^^(z) density parameter, within la and 2a errors, exist 
the possibility of being smaller compared to the corresponding noninteracting model 
(LCDM model). 

It is interesting to note that from the Figure [7| to Figure [12] and according to the 
results of the tables [H [12] and [TH the best fitted values predicted for our models are 
consistent with the requirement that the universe be older that any of its constituents 
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Figure 13. Comparison of tlie best fit reconstructed and its errors for the rate between 
dark density parameters fl*jj^{z)/n*jyj^{z) in terms of the parameters Aq, Ai, A2 and 
the corresponding Chebyshev polynomials for the models II (left above panel), III 
(right above panel) and IV (left below panel) , respectively. Red lines show best fit 
reconstructed results, while green lines and blue lines show reconstructed errors within 
the la and 21a confidence level errors, determined from a combination of SNela + 
BAO + CMB + H + X-ray dataset. In addition, note that within the la and 2a 
errors, our result and its errors are consistent with the constraints in the LCDM model. 
Furthermore, our best fit results and its errors at la and 2a are in LCDM model. 
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Figure 14. Superposition of tlie best estimates for the dimensionless interaction 
function Iq{z) (left above panel), the density parameters i}^^{z) (right above 

panel) and the deceleration parameter q{z) (left below panel) as a function of the 
redshift for the models I (green line), II (red line). III (blue line) and IV (black line). 
By comparison, the LCDM model (pink line) is shown. The curves show the best 
estimates using the expansion of all the functions in terms of the first n = 2 Chebyshev 
polynomials. Note that the reconstruction of the best estimate of the dimensionless 
interaction function Iq(z) for the models II, III and IV produces roughly the same 
curve and that the density parameter of dark energy is definite positive for all the 
range of redshift considered in the reconstruction. 



at a given redshift |102j , |92] . 

From all our analysis, and from the Figure [2] and |3l [5] and [THl show that the la 
and 2a constraints on the EOS parameter w contain high probability of being in the 
phantom region, but our results are even consistent and compatible with the LCDM 
model at lex error. 

In general when more cosmic observational data sets are added to constrain our 
model parameters space, the degeneracies between model parameters will be broken. 
The reason comes from the fact that the constraints on our models are more stringent 
that the results obtained in reference [1]. 
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7. Conclusions. 

In this paper, we developed tlieoretically a novel method for the reconstruction of the 
interaction function between dark matter and dark energy assuming an expansion of the 
general interaction term proportional to the Hubble parameter in terms of Chebyshev 
polynomials which form a complete set of orthonormal functions. To show how the 
method works, we applied it to the reconstruction of the interaction function expanding 
it in terms of only the first N Chebyshev polynomials (with = 1,2,3,4) and fitted 
for the coefficients of the expansion assuming three models: (a) a DE equation of the 
state parameter u = —1 (an interacting cosmological A), (b) a DE equation of the state 
parameter u = constant and (c) a DE equation of the state parameter u = constant 
and fi_DM = constant, respectively. The fit of the free parameters of every model is done 
using the Union2 SNe la data set from "The Supernova Cosmology Project" (SCP) 
composed by 557 type la supernovae [3J, the Baryon Acoustic Oscillation (BAO), the 
Cosmic Microwave Background (CMB) data from 7-year WMAP, the Observational 
Hubble data and the Cluster X-ray Gas Mass Fraction. 
Our principal results can be summarized as follows: 

(i) Compared with the results of reference [I], our fitting results show a faster 
convergence of the best fitted values for the several cosmological variables considered 
in this paper when the numbers of parameters A^ is increased in the expansion (H71) . 

(ii) The best estimates for the interaction function Iq(-2) prefer to cross the 
noninteracting line Iq(-2) = during the present cosmological evolution. This 
conclusion is independent of the numbers of coefficients (up to A^ = 4 in this work) 
used in the expansion of Iq{z). The crossing implies a change of sign of Iq{z) from 
positive values at the past (energy transfers from dark energy to dark matter) to 
negative values at the present (energy transfers from dark matter to dark energy). 
This direction of decay is similar to the results found in the recent literature [1] 
and is in disagreement with the oscillatory behavior reported in [88] . 




(iii) The statement above is conclusive because the existence of crossing of the 
noninteracting line Iq{z) = in some moment of the recent past is totally contained 
inside the la and 2a constraints given by current observations. 

(iv) We can state that is reasonable to expect that DE and DM interact via a small but 
calculable coupling and to be consistent with the observations at 2a error. In this 
aspect. The the Cosmic Microwave Background (CMB) data provides a stringent 
constraints on the coupling. 

(v) We confirm that adding the Baryon Acoustic Oscillation (BAO), the Cosmic 
Microwave Background (CMB), the Observational Hubble (H) and the Cluster X- 
ray Gas Mass Fraction (X-ray) data sets, generally reduces the allowed range of the 
model parameters. It is to say are strongly reduced compared with the results of 
the ref [1], which were obtained using only Union2 SNe la data. 

(vi) For all the interacting models studied in this work, the best estimates for the dark 




Dark energy interacting with dark matter 



33 



energy density parameter Q*jj^{z) becomes positive definite in the range of redsfiifts 
considered in tliis work. Tliis statement is conclusive because, witliin tfie la and 
2a errors for the fit with three parameters (A^ = 2), Q*j^^{z) becomes positive in all 
the range of redshifts considered in the observational data sets. 

(vii) The la and 2a confidence intervals, for the EOS parameter u considered in the 
marginalized models III and IV, preferred to be in the phantom region {w < —1). 
But exists a small probability for the EOS parameter u of being in the quintessence 
region. This conclusion is consistent and is within of LCDM model. 

(viii) According to the Figures [6] and [11] and to our results, we confirm that there is strong 
evidence that the universe is accelerated in the recent past, which is consistent with 
studies in [I25], [T26] . 

(ix) we therefore conclude that our models are consistent with the properties of the 
cosmological variables and its results are also consistent with the LCDM model at 
2a error. 

Finally, we emphasize the importance of more accurate and extensive 
measurements. These observations will certainly provide a complementary tool to test 
the reality of the current cosmological parameters, of the term of interaction between 
dark sectors and the existence of the crossing of the noninteracting Iq(;z) = line and 
place stringent constraints on them, and as a consequence, to distinguish among the 
many alternatives models. We believe that the analysis, combined with the results give 
us the possibility of understanding all that. With the successful experience achieved 
to reconstruct the general interaction term between dark matter and dark energy, our 
next step is to extend our study to do a double reconstruction as much the interaction 
term Iq{z) as dark energy equation of state parameter u simultaneously, within the dark 
sectors. More efforts will be required on this way to carry out it. 
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